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Lateral heterostructures of two-dimensional materials may exhibit various intriguing emergent 
properties. Yet when specified to the orientationally aligned heterojunctions of zigzag graphene and 
hexagonal boron nitride (/iBN) nanoribbons, realizations of the high expectations on their properties 
encounter two standing hurtles. First, the rapid accumulation of strain energy prevents large- 
scale fabrication. Secondly, the pronounced half-metallicity predicted for freestanding graphene 
nanoribbons is severely suppressed. By properly tailoring orientational misalignment between zigzag 
graphene and chiral /iBN nanoribbons, here we present a facile approach to overcome both obstacles. 
Our first-principles calculations show that the strain energy accumulation in such heteroj unctions 
is significantly diminished for a range of misalignments. More strikingly, the half-met allicity is 
substantially enhanced from the orientationally aligned case, back to be comparable in magnitude 
with the freestanding case. The restored half-met allicity is largely attributed to the recovered 
superexchange interaction between the opposite heteroj unction interfaces. The present hndings 
may have important implications in eventual realization of graphene-based spintronics. 

PACS numbers: 81.10.Aj, 73.20.-r, 73.22.-f, 61.46.-w 


Recent research in the design, synthesis, and property 
characterization of two-dimensional (2D) heterostruc¬ 
tures represents a major advance in low-dimensional ma¬ 
terials science [T]. As compelling examples, lateral het¬ 
erostructures of graphene (G), hexagonal boron nitride 
(/iBN), and other 2D materials [21JT4] have been suc¬ 
cessfully fabricated, where the 2D materials with distinc¬ 
tively different band structures are integrated by lateral 
covalent bonding within a single atomic layer. This new 
class of hybrid materials may exhibit novel and diverse 
properties [Mil], which are expected to have a broad 
range of applications in electronic devices. 

To date, extensive experimental and theoreti¬ 

cal dZlllg] investigations have been devoted to orienta¬ 
tionally aligned G-hBN heteroj unctions. Whereas such 
heterostructures have clear advantages, there exist stand¬ 
ing challenges in actual realization of the intriguing elec¬ 
tronic and spintronic properties of graphene nanoribbons 
(GNRs). The first is related to the fabrication of the G- 
/iBN heterostructures. It has been shown experimentally 
that small pieces of orientationally aligned G-hBN lateral 
heteroj actions can be formed with high degrees of crys¬ 
tallinity and abrupt interfaces isHini, but such coherently 
strained samples cannot be enlarged significantly, mainly 
due to their lattice mismatch. It is known that in epi¬ 
taxial growth, the strain caused by lattice mismatch is 
nonlocal m, and the associated strain energy accumu¬ 
lates along with the domain size. As such, the release of 
strain energy upon a critical domain size proceeds via ini¬ 
tiation of structural instabilities such as dislocations mi. 


The other major and more fundamental challenge is as¬ 
sociated with the preservation of the electronic and spin¬ 
tronic properties. Specifically, freestanding zigzag GNRs 
have been predicted to exhibit half-metallicity under a 
strong transverse electric field [21] |22] . Because it has not 
been possible to apply the required strong external fields 
at such nanoscales, alternative approaches have been pro¬ 
posed to establish the required fields internally, for exam¬ 
ple via proper molecular adsorption [23] or sandwiching 
a GNR between two hBN NRs p^SHTTl ITO] . However, 
even though the electric field associated with the charge 
transfer between the two inequivalent G-hBN boundaries 
is indeed strong mi, the corresponding half-metallicity 
has been shown to be severely suppressed from that of 
freestanding GNRs [24]. 

In this Letter, we provide a new scheme to overcome 
the main challenges described above, by invoking orien¬ 
tationally misaligned lateral heteroj unctions consisting 
of zigzag G and chiral hBN NRs that otherwise satisfy 
structural commensurability conditions [25]-[28]. Using 
first-principles calculations, we show that the strain en¬ 
ergy accumulation in such heteroj unctions can be dras¬ 
tically diminished regardless of the domain size. More 
importantly, we reveal substantial enhancements in the 
half-metallicity from the orientationally aligned case, 
back to be comparable in magnitude with that of the 
freestanding GNR. Through detailed analysis of the spa¬ 
tial charge distribution, we attribute the restored half- 
metallicity to the enhancerecoveredd superexchange in¬ 
teraction, which greatly reinforces the coupling between 
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FIG. 1 (color online). Lateral heterojunctions of zigzag 
G and /iBN NRs with (a) orientational alignment and (b) 
orientational misalignment, where the GNRs contain iV = 4 
zigzag G chains. The balls represent the atoms inside the 
computational supercell. For each G chain, there are m = 2 
and 20 G atoms in the unit cell for the (a) orientational aligned 
and (b) misaligned case, respectively. Lattice vectors ai and 
a 2 (bi and b 2 ) define the hexagonal lattice for /iBN (G). In 
(a) the dashed lines guide the unit cell. In (b) the dislocations 
formed at the interfaces are grey shaded, (c) Accumulations 
of the strain energy as a function of the GNR width measured 
by the number of zigzag G chains, N, for both orientationally 
aligned and orientationally misaligned hetero junctions. 

TABLE 1. Possible solutions to Eq. 0 under the constrain 
of T 2 = (I 5 O) for lateral heterojactions consisting of zigzag 
GNR and chiral hBN NR. _ 

h I2 (ni,n2) 6 {%) on ~ 


1 

10 

(3,8) 

0.26 

15.3 

1 

11 

(3,9) 

0.10 

13.9 

1 

12 

(3, 10 ) 

0.02 

12.7 

1 

13 

(3,11) 

0.02 

11.7 

1 

14 

(3,12) 

0.03 

10.9 

1 

15 

(3,13) 

0.03 

10.2 


the spin states at the two opposite interfaces. 

The concept of commensurability matching has been 
widely invoked in studies of grain boundaries [25j |26l [28] 
and heterostructures [27] in 2D materials. For the zip¬ 
ping of zigzag G and chiral /iBN NRs, it simply reads, 

/iTi=/2T2, (1) 

where, Ti and T 2 are the primitive translation vectors 
along the interface for chiral /iBN NR and zigzag GNR, 
respectively, as shown in Figs. 1(a) and 1(b). Here li and 
I 2 are integer numbers, with li = 1 for most cases. Using 
the lattice vectors shown in Fig. 1(b), it is convenient to 


express Ti = niai + n 2 a 2 or (ni,n 2 ) and T 2 = bi or 
(1,0), with |Ti| = + nin2 + n 2 ^ and IT2I = bo, 

ni and n 2 are integer numbers, while ao = 1.450 A and 
bo = 1.425 A are the length of the lattice vectors for /iBN 
and G, respectively. The orientational misalignment be¬ 
tween /iBN and G also equals the chiral angle of /iBN 
NRs, which is defined as, 

0 = tan“^ [\/ 3 n 2 /( 2 ni + 77.2)]. ( 2 ) 

For the specific case shown in Fig. 1(b), Ti = (3,8), 
T 2 = (1,0), and 0 = 15.3°. In general, Eq. 0 can only 
be approximately solved such that a residue value in the 
lattice mismatch still remains, 

<^= l^ilTil-;2|T2||/ai|Ti|). (3) 

Table lists the optimal solutions to Eq. 0 that mini¬ 
mize 6 for a series of possible misalignment angles. We 
note that for the orientationally aligned case as shown in 
Fig. 1(a), the lattice mismatch is relatively large, given 
hy 6 = (ao — bo)/ao ^ 1.8%. In contrast, here it can 
be seen that the lattice mismatch is much smaller, with 
5 < 0.3%. The vanishingly small lattice mismatch indi¬ 
cates that the commensurability conditions are well sat¬ 
isfied. We therefore can expect that the accumulation of 
strain energy in the orientationally aligned heterostruc¬ 
ture can now be essentially avoided. This is of critical 
importance for the fabrication of large-scale G-ZiBN het¬ 
erostructures that may exhibit various intriguing elec- 
tronic/spintronic properties. 

To quantitatively confirm these expectations on the en¬ 
ergetic and electronic properties, we have carried out sys¬ 
tematic first-principles density functional theory (DFT) 
studies of various G-ZiBN heterojunctions, including both 
orientationally aligned and misaligned cases. For the 
misalignment cases, we choose the heterostructures with 
the same misalignment angle of 0 ~ 15.3° as shown in 
Fig. 1(b), but with varying GNR widths. We also note 
that this case corresponds to the largest residue mis¬ 
match, with all the other misalignments given in Table [l| 
exhibiting substantially lower residue mismatches. Our 
DFT calculations were made using the Vienna ab initio 
simulation package (VASP) [29], where the projector aug¬ 
mented plane wave (PAW) method [30l|3T] was used, and 
the generalized gradient approximation (GGA) in the 
framework of Perdew-Burke-Ernzerhof (PBE) [32] was 
adopted for the exchange-correlation interaction. The 
positions of the atoms were obtained by structural opti¬ 
mization until the forces on each atom are smaller than 
0.01 eV/A. We have also taken the corrugations associ¬ 
ated with out-of-plane atomic displacements into consid¬ 
eration, and used a vacuum layer of 20 A to avoid the 
possible effects of image supercells. A plane-wave basis 
was set with a kinetic-energy cutoff of 500 eV, and the 
Brillouin zone was sampled bya5x5xlor5x31xl 
/c-mesh, depending on the supercell size. 
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FIG. 2 (color online). The spin-resolved band structures of 
orientationally (a) aligned and (b) misaligned lateral hetero¬ 
junctions of zigzag G and /iBN NRs for GNRs with different 
widths, N. The T ^ X direction corresponds to the inter¬ 
facial direction spatially. The misalignment angle in (b) is 
0 ^ 15.3°. The solid (blue) and dashed (red) curves represent 
different spin orientations. In (b) the half-metallic gap is grey 
shaded. The Fermi energy is set at zero as indicated by the 
horizontal short dashed lines. 

For an effective characterization of the strain energy 
accumulation in the GNR, the length of the supercell par¬ 
allel to the interfacial direction is determined by the lat¬ 
tice constant of pristine /iBN and fixed during the struc¬ 
tural optimization, while the length perpendicular to the 
interfacial direction is optimized energetically. We calcu¬ 
late the strain energy versus the GNR width in terms of 
the number of zigzag C chains, A^, as, 

E = Ejsf — F/2 — tti[N — 2 )Eq^ (4) 

where or E 2 represents the formation energy of the 
heterostructure with N or 2 carbon chains, respectively, 
Eg denotes the formation energy of pristine graphene 
per C atom, and m is the number of C atoms in one C 
chain contained in the supercell. Note that m = 2 and 
20 for the heterostructures shown in Figs. 1(a) and 1(b), 
respectively, while Ejy, E 2 , and Eq are all obtained from 
detailed first-principles calculations. In particular, the 
term {E^ — E 2 ) cancels out not only the contributions 
from the respective /iBN components in orientationally 
aligned and misaligned cases, but also the contributions 
from the formation of the C-B and C-N bonds at both 
interfaces for either alignment. This choice ensures that 
the strain energy E contains purely the energy raised 
from the enlargement or “growth” of the GNR width. 
Finally, E is normalized to energy per unit length along 
the interfacial direction. 

The results are summarized in Fig. 1(c). First of all, 
for both the orientationally aligned and misaligned cases, 
the strain energy increases rapidly for N < 4, largely 


TABLE II. Gomparison of the half-metallic gaps (meV) of ori¬ 
entationally aligned heterojunctions, misaligned heterojunc¬ 
tions, and freestanding GNRs at different GNR widths, N. 
The results in Ref. m and Ref. [22] are shown as references. 


N 

Aligned 


Misaligned 

Ereestanding GNR 


This work Ref. [17] 

PBE PBE 

This work 

PBE 

i’his work 
PBE LDA 

Ref. Ea 

LDA 

8 

75 

71 

160 

691 

405 

443 

10 

18 

0 

160 

629 

362 

- 

12 

0 

0 

138 

555 

327 

- 

14 

0 

- 

131 

471 

266 

- 

16 

0 

- 

112 

425 

276 

290 

18 

0 

- 

106 

386 

271 

- 


due to the local structural relaxation near the interfaces. 
In the N > 4 region, distinctly different behaviors of 
the strain energy are identified. For the orientationally 
aligned case, the strain energy increases monotonously 
with a near-linear dependence on N, indicating rapid 
accumulation of the strain energy. This finding is con¬ 
sistent with the recent experimental revelation m that 
defects such as dislocations, discontinuities, and ripples 
emerge to release the accumulated strain energy in the 
heterostructure of G-hBN during the epitaxial growth. 
On the other hand, for the orientationally misaligned 
case, it is reassuring to see that the strain energy in¬ 
deed stays at a relatively low constant level, independent 
ofN. 

The preceding analysis reveals that the superior struc¬ 
tural stability can be achieved in properly misaligned het¬ 
erojunctions of zigzag G and chiral /iBN NRs, where the 
accumulation of strain energy can be effectively dimin¬ 
ished. As a consequence, hetero junctions at relatively 
large scales can be fabricated. The next question is how 
such intentional misalignments would influence the phys¬ 
ical properties of the heterostructures, which contain pe¬ 
riodic arrays of dislocations in the form of pentagon- 
heptagon pairs along the interfaces [see Fig. 1(b)]. In 
this regard, we note that similar structural topologies 
in polycrystalline graphene have been found to enhance 
the mechanical stability [33| , to alter the transport prop¬ 
erties across the grain boundaries [25|, or to introduce 
emergent magnetic properties along the boundaries [34] . 

Figure 2 compares the band structures of orientation¬ 
ally aligned and misaligned G-hBN lateral heterojunc¬ 
tions for GNRs with different widths. Although the 
band structures are spin-polarized, and half-metallicity 
can be identified for both cases, the differences between 
the two are substantial. For the aligned cases, the band 
structures display half-metallicity only for N < 10. The 
half-metallic gap decays rapidly with and is already 
nearly closed when N > 10, indicating an effective con¬ 
ventional metal, Fig. 2(a). On the other hand, for the 
misaligned cases, the band structures reveal pronounced 
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FIG. 3 (color online). Spatial charge distributions of the edge 
electronic states within the energy range of 0.0 — 0.5 eV below 
the Fermi level, for both orientationally (a) aligned and (b) 
misaligned lateral heterojunctions of zigzag G and /iBN NRs. 
The isosurface level is 0.01 e/A . The misalignment angle in 
(b) is ^ ~ 15.3°. (c) and (d) are the side view of the (a) and 
(b), respectively, displaying the vertical corrugations in the 
misaligned case. The black dashed lines guide the interfaces. 


half-metallicity for all the range of N considered. The 
half-metallic gap is relatively large, and the decay with 
N is not as significant, Fig. 2(b). We have also calcu¬ 
lated quantitatively the half-metallic gaps of freestanding 
zigzag GNRs. The results are summarized in Table |TI[ 
Compared to freestanding GNRs, the half-metallic gap is 
severely suppressed in the orientationally aligned hetero¬ 
junctions, while restored greatly in the misaligned het- 
erojunctions. These results can be rationalized by com¬ 
paring with previous reports [nnn]. For orientationally 
aligned hetero junctions and freestanding GNRs, our re¬ 
sults show good agreement with those from the earlier 
studies. Of course, it should be noted that the PBE re¬ 
sults are ^ 40% larger than the local density approxima¬ 
tion (LDA) results, as the PBE functional overestimates 
the band gaps, while the LDA functional underestimates 
the band gaps. 

The enhanced half-metallicity in the misaligned cases 
originates from the enhanced superexchange interaction 
between the opposite interfaces. This insight is obtained 
by analyzing the spatial charge distribution along the in¬ 
terfaces. Recall that in freestanding zigzag GNR, the 
edge states are spin-polarized with ferromagnetic order 
along a given edge but are antiparallel for the two oppo¬ 
site edges |35ti38] . These features are critical for achiev¬ 
ing half-metallicity. In particular, due to the localization 
of the electronic density at either side of a freestand¬ 
ing GNR (35] |36] , the interedge superexchange interac¬ 
tion is pronounced [3S1I10]. This superexchange interac¬ 
tion opens a gap for the ground states, and at the pres¬ 
ence of a strong external electric field, a sizeable half- 
metallic gap can be generated. Unfortunately, such ideal 
spatial charge distributions cannot be readily preserved 
for the heterojunctioned systems. Eor the orient at ion- 


TABLE III. Gomparison of AE (meV/nm) favoring AFM 
in the orientationally aligned heterojunctions, misaligned 
hetero junctions, and freestanding GNRs at different GNR 
widths, N. 


N 

Aligned 

Misaligned 

Freestanding GNR 

8 

4 

10 

335 

10 

2 

8 

341 

12 

1 

7 

352 

14 

1 

6 

334 

16 

1 

6 

367 

18 

1 

6 

372 


ally aligned cases, the charge density of the edge states 
at one interface is vanishingly small, Eig. 3(a), presum¬ 
ably due to the asymmetric bonding of C-B and C-N 
at the opposite interfaces [24]. Consequently, the su¬ 
perexchange interaction between the opposite interfaces 
is drastically diminished. In contrast, for the orienta¬ 
tionally misaligned cases, the charge of the edge states 
is distributed inhomogeneously along the two interfaces 
due to the presence of the dislocations, Eig. 3(b). The 
charge densities in the misaligned cases are lower than 
that of freestanding GNRs [351438] . but are comparable 
in magnitude along the two opposite interfaces, leading 
to the largely restored strengths in the superexchange 
interaction and half-metallicity. 

The enhanced superexchange interaction in the orien¬ 
tationally misaligned cases is further supported by the 
energy difference between the anti-ferromagnetic (AEM) 
and non-magnetic (NM) state, 

AE = Enm - ^AFM, (5) 

where Eafm and Enm denote the energy of the AEM and 
NM states, respectively. The calculated AE is scaled to 
energy per nanometer in the interfacial direction. Larger 
AE corresponds to stronger superexchange interaction. 
Thus, according to the data shown in Table |III[ the 
superexchange interaction is strongest for freestanding 
GNRs, and is stronger for the orientationally misaligned 
cases than the orientationally aligned cases. 

We now briefly discuss the possibilities of experimen¬ 
tally realizing orientationally misaligned G-/iBN lateral 
heteroj unctions. We notice that when graphene is grown 
first, the /iBN epilayer will take the crystallographic ori¬ 
entation of the G seed layer uni. On the other hand, it is 
feasible to obtain orientational misalignment between G 
and /iBN by exploiting the different adhesive properties 
between G and /iBN to a specific catalytic substrate such 
as Cu(lOO). It has been shown that the growth of hBN 
undertakes definitive crystallographic orientations follow¬ 
ing the square lattice of the Cu(lOO) surface [41], while 
the growth of G assumes rather random crystallographic 
orientations [lol EH ES]. Therefore, a specific orienta¬ 
tional misalignment between G and /iBN can be achieved 
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by first selecting G with a proper crystallographic orien¬ 
tation from the randomly distributed samples grown on a 
Cu(lOO) substrate. Next, hBN can be grown at different 
areas on the Cu(lOO) substrate using mask approaches, 
and eventually coalesces with the G domain forming a 
misaligned G-hBN junction. Overall, given the clear ad¬ 
vantages of such misaligned heterojunctions, we believe 
other more creative approaches will be devised on the 
experimental side for their realization. 

In summary, we have introduced a new way to form 
zigzag G-hBN hetero junctions by properly tailoring the 
orientational misalignment between the two. The strain 
energy accumulation in such misaligned heterostructures 
is essentially eliminated, while the half-metallicity is 
found to be drastically enhanced from the orientation- 
ally aligned cases, back to be comparable in magnitude 
to that of freestanding GNRs. The pronounced half- 
metallicity is further attributed to the restored strength 
of the superexchange interaction between the electrons 
located at the two opposite interfaces. Overall, these 
revelations are valuable for potential practical realiza¬ 
tion of the intriguing emergent electronic and spintronic 
properties of G-hBN heterostructures. 
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